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Abstract: For a small library of natural products from marine sponges and ascidians, 
in silico docking to the Lymnaea stagnalis acetylcholine-binding protein (AChBP), a 
model for the ligand-binding domains of nicotinic acetylcholine receptors (nAChRs), was 
carried out and the possibility of complex formation was revealed. It was further 
experimentally confirmed via competition with radioiodinated a-bungarotoxin ([ I]-aBgt) 
for binding to AChBP of the majority of analyzed compounds. Alkaloids pibocin, varacin 
and makaluvamines C and G had relatively high affinities {Ki 0.5-1.3 |a,M). With the 
muscle-type nAChR fi-om Torpedo californica ray and human neuronal a7 nAChR, 
heterologously expressed in the GH4C1 cell line, no competition with [^^^I]-aBgt was 
detected in four compounds, while the rest showed an inhibition. Makaluvamines 
{Ki ~ 1.5 |a,M) were the most active compounds, but only makaluvamine G and 
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crambescidine 359 revealed a weak selectivity towards muscle-type nAChR. Rhizochalin, 
aglycone of rhizochalin, pibocin, makaluvamine G, monanchocidin, crambescidine 359 
and aaptamine showed inhibitory activities in electrophysiology experiments on the mouse 
muscle and human a7 nAChRs, expressed in Xenopus laevis oocytes. Thus, our results 
confirm the utility of the modeling studies on AChBPs in a search for natural compounds 
with cholinergic activity and demonstrate the presence of the latter in the analyzed marine 
biological sources. 

Keywords: marine natural compounds; acetylcholine-binding protein; nicotinic 
acetylcholine receptors; computer modeling; radioligand assay; electrophysiology 



1. Introduction 

Nicotinic acetylcholine receptors (nAChRs) are among the most comprehensively studied 
ligand-gated ion channels (see reviews [1-4]). They are widespread in fish electric organs and 
mammalian muscles (so-called muscle type of nAChR), central and peripheral nervous systems 
(neuronal nAChRs) and some non-neuronal tissues (so-called "non-neuronal" nAChRs), where they 
perform a multitude of functions fi-om conducting nerve-muscle transmission to participation in the 
different cognitive processes and regulation of inflammatory response. Involvement of distinct nAChR 
subtypes in different pathologies (muscle dystrophies, Alzheimer's and Parkinson's diseases, 
schizophrenia, nicotinic addiction, chronic pain) dictates the need for potent and highly selective 
cholinergic ligands to the respective receptor subtypes. In addition, utilization of similar new 
compounds as drugs in clinics is under intensive development [5]. 

A huge number of compounds of different chemical nature fi-om various taxa (bacteria, plants, 
moUusks, chordates) are known to interact with nAChRs. Among them, the most well-known and 
widely used are a-neurotoxins from snake venoms (see reviews [6-8]). The last three decades 
witnessed the discovery and intensive scientific application of numerous Conus moUusk 
venom-derived conotoxins of various classes (see reviews [9-11]). Other marine creatures were much 
less studied for presence of cholinergic compounds remaining in the shadow of the latter. Among 
them, spirolide and gymnodimine phytoplankton toxins [12], nereistoxin from Lumbriconereis heteropoda 
annelid [13] and two ascidian alkaloids [14] should be mentioned. 

Herein, we describe cholinergic properties of 13 natural low molecular weight compounds, isolated 
at PIBOC from marine sponges and ascidians (Figure 1). A large number of various bioactive 
compounds were earlier isolated from these two animal taxa, however molecular targets were not 
identified for the most of them. For some of them, a structural similarity to diverse cholinergic ligands 
(quaternary ammonium salts, heterocyclic compounds) allowed us to anticipate their possible activity 
towards nAChRs. To check this, we performed docking of these natural products to Lymnaea stagnalis 
acetylcholine-binding protein (AChBP) using the available X-ray structures of this protein in 
complexes with different cholinergic ligands. Several known AChBPs were found to be excellent 
structural models for the ligand-binding domains of all nAChRs (see reviews [15,16]), and now are 
widely used from purification of new natural cholinergic ligands [17] to design the libraries of 
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synthetic compounds [18,19]. In the present communication, the conclusions from computer modeling 
were verified by efficient interaction of the studied compounds with AChBP revealed by radioligand 
analysis, as well as by their binding to muscle and al neuronal nAChRs tested by radioligand analysis 
and elecfrophysiology. 

Figure 1. Chemical structures of compounds from marine sponges and ascidians (1-13), 
for which putative cholinergic activities were examined by computational and 
experimental methods. 




2. Results and Discussion 

2.1. Isolation of Individual Compounds 

Structures of the tested compounds are given on the Figure 1 . 

Rhizochalin (1), its aglycone (2) [20,21], pibocin (3) [22], and monanchocidin (13) [23] were for 
the first time isolated and described as new natural products by authors of the paper (T.M. and 
colleagues) from the sponge Rhizochalina incrustata, an ascdian Eudistoma sp. and the 
sponge Monanchora pulchra, respectively. Varacin (4) was isolated from the far-eastern ascidian 
Polycitor sp. [24]. Makaluvamines C and G (5, 6) were isolated from the Ausfralian sponge 
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Zyzzya fuliginosa [25] and structurally identified with the corresponding alkaloids previously obtained 
from Fijian sponges belonging to the same genus [26], and an hidonesian sponge Histodermella sp. [27]. 
Debromohjmienialdesine (7) was isolated from the sponge Acanthella cartery [28] and identified with 
the alkaloid of the sponge Phakellia flabelata [29]. l,r-Dimethyl-[2,2']-bipyridyldiium salt (8) was 
isolated from the Far-Eastern ascidian Botrylloides violaceus and identified with the same compound 
found in the bivalve moUusk Callista chione earlier [30]. Crambescidin 359 (9) was obtained from 
alcoholic exfract of the Ausfralian sponge Monanchora clathrata and identified with the alkaloid 
previously isolated from the sponge Monanchora unguiculata [31]. 7,8-Dihydroimidazo-[l,5-c]- 
pyrimidin-5(6//)-one (10) was isolated from the Vietnamese sponge Aplysina sp. [32] and structurally 
identified with the metabolite of Aplysina fistularis forma fulva [33]. 1,3-Dimethylisoguaniniium 
hydrochloride (11) was isolated from the Far-Eastern ascidian Syncarpa oviformis and identified with 
the same compound from the sponge Amphimedon paraviridis [34]. Aaptamine (12) was obtained 
from extracts of the Vietnamese sponge Aplysina sp. [35] and identified at comparison with that from 
Aaptos aaptos [36]. 

2.2. Docking toAChBP and Analysis of Binding Parameters in Competition with [^^^IJ-aBgt 

We performed docking of all compounds to Lymnaea stagnalis AChBP using the structure of a 
protein with a HEPES buffer molecule in the binding site [37], the AChBP structure in complex with 
competitive antagonist dihydro-P-erythroidine (DHpE) [38], which is a low molecular weight alkaloid, 
and with the competitive antagonist a-cobratoxin [39], a snake a-neuro toxin close in structure and 
properties to a-bungarotoxin. In the a-cobratoxin complex, the loop C of each protomer moved away 
from the central axis of pentameric AChBP in comparison with its position in the HEPES bound protein. 
On the contrary, in the case of the complex with DH^E, its position is even more close to the central axis 
and is shifted down in the perpendicular direction. In our opinion, the use of three different spatial forms 
of AChBP for docking each of the studied compounds increases the significance of such calculations. 

The results, presented in the form of calculated inhibition constants {Ki), are summarized in Table 1 . 
For almost all the compounds, the most energy-preferable position was found in the "classical" binding 
site corresponding to the hydrophobic pocket of nAChRs, where the cholinergic agonists and 
competitive antagonists are bound. However, computer docking in some cases resulted in the preferred 
location of the molecule in other (non-classical) binding sites. The most sfriking example is the 
1,3-dimethylisoguaniniium hydrochloride (11), for which the theoretical calculated Ki values for 
"classical" and another binding sites differed by a factor of 8 (see Table 1). 

According to computational data, a high affinity to L. stagnalis AChBP could be expected for some 
compounds (3-7, 9, 12, 13), and low efficiency of interaction or lack of it — for the others 
(compounds 1, 2, 8, 10, 11). It should be noted that docking of almost all active compounds (3, 4, 5, 6, 
9, 13) to three different forms of AChBP gave similar results with predicted affinities differing by less 
than 5 |j,M. Only in case of the moderately active aaptamine (12) and very weak l,r-dimethyl-[2,2']- 
bipyridyldiium salt (8) the predicted affinities displayed significant discrepancy, which could possibly 
be explained by differences in distances between the AChBP residues critical to a-cobratoxin binding, 
as maximal differences are observed in docking to that form of AChBP (Table I, column C). 
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Table 1. The calculated and measured affinities of the studied compounds (see Figure 1) to 
L. stagnalis acetylcholine-binding protein (AChBP). A — ^HEPES-bound; B — DHpE-bound 
and C — a-cobratoxin-bound forms of AChBP respectively; * — only a part of the 
compound could be docked correctly because of multiple rotatable bonds; ** — ^the binding 
outside "classical" agonist/competitive antagonist site is predicted; n.d. — ^not determined (no 
positive solutions). The measured affinities were obtained as follows: firstly, the respective 
IC50 values were calculated from inhibition curves (see Figure 2) using ORIGIN 7.5 and 
then converted into the X,s according to Cheng-Prusoff equation (see Experimental 
Section) with the mean ± s.e.m. values of duplicate or triplicate measurements for each 
compound concentration. 

Theoretical calculated Ki values (jiM) for 

Experimental Ki values (jiM) at interaction 

Compounds interaction with L. stagnalis AChBP 

with L. stagnalis AChBP from radioligand assa 



1 * 


n.d. 


n.d. 


39 


28 ±4 


2 * 


n.d. 


n.d. 


39 


130 ± 10 


3 


0.98 


1.7 
0.64 ** 


2.0 


0.83 ± 0.04 


4 


1.6 


0.97 


0.75 


0.79 ± 0.03 


5 


2.2 


0.46 


5.8 


1.3 ±0.1 


6 


0.05 


0.04 


0.97 


0.55 ±0.01 


7 


8.2 


0.79 
0.50 ** 


1.7 


>1000 


8 


33 


12 


310 


540 ± 60 


9 


0.08 


0.47 


0.40 


27 ±2 


10 


170 


120 


n.d. 

320 ** 


>1000 


11 


150 


83 
J J ** 


350 


>1000 


12 


1.5 


5.3 


19 


3.0 ±0.4 


13 * 


0.53 


0.43 


0.11 


2.5 ±0.4 



The interactions, suggested by the calculations, were confirmed in radioligand assay by testing the 
ability of the compounds to compete with [^^^I]-aBgt for the L. stagnalis AChBP "classical" binding 
sites. The inhibition curves of the most active compounds are presented in Figure 2 and the respective 
Ki values are collected in Table 1 . The compounds with the worst in silico docking score showed also a 

125 

poor ability to compete with [ I]-aBgt (1, 2, 8) or were inactive with Ki > 1000 [iM (compounds 10, 11). 
Moreover, binding activities for some of the most active compounds (3-6, 12, 13) showed good 
agreement with the results of computer docking. 

Thus, the results obtained for the majority of the studied compounds show once again sufficient 
validity of this computational approach for a preliminary evaluation of affinity of new ligands and 
even for generating assumptions about intermolecular bonds determining this affinity. 

As an example, we can consider a model of the L. stagnalis AChBP complexed with the compound 
that shows the best binding parameters both in the computer docking and in radioligand assay, namely 
with makaluvamine G (6) (Figure 3 a). In this schematic planar model all the ligand interactions with 
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formation of hydrogen bonds and hydrophobic contacts with detinite AChBP amino acid residues can 
be clearly seen. A high affinity of makaluvamine G (6), deduced from computer calculations, is 
determined by two hydrogen-bonds and 14 hydrophobic contacts with different L. stagnalis AChBP 
amino acid residues (Figure 3 a). 

125 

Figure 2. Inhibition of [ I]-aBgt binding to L. stagnalis AChBP with the most active 
compounds studied. Numbering the compounds corresponds to the numbering in Figure 1. 
The corresponding curves are marked with symbols: 1 — filled circles; 2 — open circles; 
3 — filled squares; 4 — open squares; 5 — filled diamonds; 6 — open diamonds; 
9 — open triangles; 12 — filled triangles; 13 — stars. Each point is a mean ± s.e.m value of 
two or three measurements for each concentration. The curves were calculated from the 
means ± s.e.m. using ORIGIN 7.5 program (see Experimental Section). The respective Kj 
values are listed in Table 1 . 

100- 
Cjb75- 

■3 

•9 50- 
w 

•y 25- 
V 

a 



8 7 6 5 4 3 2 

- Ig [compound, M] 

We compared this calculated structure with the X-ray structure of the same protein complexed with 
clothianidin (a neonicotinoid insecticide) (Figure 3b) [40]. In that case, two hydrogen bonds (one of 
which with the oxygen of bound water molecule) and 10 hydrophobic interactions were revealed, 6 of 
which are in contacts with the same AChBP amino acid residues as in the complex with makaluvamine G. 
The similarity of the location of these two molecules in the binding site is obvious. Note that the 
affinity of clothianidin to L. stagnalis AChBP {Kd 13 |j,M evaluated by quenching of intrinsic 
tryptophan fluorescence) [40] was lower than that of makaluvamine G (6) {Ki = 0.55 ± 0.01 |J.M), 
perhaps because makaluvamine G fits binding site more tightly, forming more of hydrophobic contacts 
(Figure 3c). 

Nevertheless, in two cases (compounds 7, 9) an essential discrepancy between the calculations and 
experimental data was found (see Table 1). Of course, using computer methods of calculation, one 
should never forget about false-positive solutions. However, another explanation of the disagreement 
between the calculated and measured affinities might be the interactions of tested compounds outside 
the "classical" agonist/competitive antagonist binding site. Such a possibility was simulated for the 
compound (7) (see Table 1). By comparing the in silico docked debromohjmienialdesine (7) and the 
structure of three-finger toxin a-cobratoxin in their complexes with AChBP (Figure 3d) we could 
expect poor competition of this compound with a-bungarotoxin (very similar to a-cobratoxin), because 
of partially overlapping their binding sites on this target. In fact, [ I]-aBgt binding to L. stagnalis 
AChBP could be observed in radioligand assay even in the presence of debromohjmienialdesine (7) as 




Mar. Drugs 2014, 12 



1865 



predicted by docking simulations. Analogously, an interaction of compound (9) outside the "classical" 
binding site could be supposed because of its functional activity towards al nAChR subtype (see 
Section 2.4.) which is the most close in pharmacological profile to L. stagnalis AChBP. 

Figure 3. Schematic planar image of the model of makaluvamine G (6) docked to 
HEPES-bound form of L. stagnalis AChBP binding site (a) in comparison with the same 
presentation of the X-ray structure of the complex of the same protein with clothianidin 
(PDB ID — 2ZJV) (b); intermolecular bonds in ligands are colored in magenta; the same 
links in the AChBP amino acid residues involved in the formation of hydrogen bonds 
{green lines) are shown in orange. The AChBP amino acid residues forming hydrophobic 
contacts with ligands are presented as red combs. Those of them that are the same for both 
ligands are encircled by red lines. Atoms of carbon, oxygen, nitrogen, and sulfur are 
colored in black, red, blue and yellow, respectively. Water molecule involved in the 
formation of hydrogen bonds is shown in cyan; (c) superposition of 3D structures of 
makaluvamine G (6), docked to AChBP and crystal structure of clothianidine AChBP 
complex. Carbons of makaluvamine G (6), clothianidine and AChBP are shown in light 
blue, pink and grey respectively; oxygens are shown red, nitrogens are shown blue; 
hydrogen bonds are brown; (d) superposition of a-cobratoxin structure (blue) 
(PDB ID — 1YI5) and the best solution for debromohymenialdesine (7) (magenta) in silico 
docked to L. stagnalis AChBP (gray). 




2.3. Analysis of Competition with [ IJ-aBgt for Interactions with the Torpedo Californica and a7 nAChRs 

The same radioligand assay, based on competition with [^^^I]-aBgt, was used to evaluate the affinity 
of tested compounds towards muscle-type nAChR fi-om T. californica ray and to human al nAChR 
heterologously expressed in the GH4C1 cell line. Inhibition curves for the most potent compounds and 
their Ki values are presented in Figure 4 and Table 2, respectively. 
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Figure 4. Inhibition of initial rate for [ I]-aBgt binding to T. californica nAChR (a) or 
human a7 nAChR (b) with the most active compounds studied. Numbering the compounds 
and their symbols correspond to the numbering in Figure 1 and symbols in Figure 2, 
respectively. Each point is a mean ± s.e.m value of two or three measurements for each 
concentration. The curves were calculated from the means ± s.e.m. using the ORIGIN 7.5 
program (see Experimental Section). The respective Ki values are listed in Table 2. 



(a) (b) 




765437 6 S 4 3 

- ]g [ compoimd,M] - Ig [compound, M] 



Table 2. Affinity of all the studied compounds (see Figure 1) tested in competition with 

125 

[ I]-aBgt for binding to muscle-type nAChR from T. californica and human neuronal a7 
nAChR (see respective inhibition curves in Figure 4). The IC50 values were calculated 
using ORIGIN 7.5 and then converted into the KiS according to Cheng-Prusoff equation 
(see Experimental Section) with the mean ± s.e.m. values of duplicate or triplicate 
measurements for each compound concentration. 



Compounds 


^/Values (fiM) on... 


T. californica nAChR 


human a7 nAChR 


1 


22 ±2 


33 ±2 


2 


24 ±2 


56 ±4 


3 


44 ±7 


130 ±8 


4 


10 ± 1 


19 ± 1 


5 


1.30 ±0.05 


3.6 ±0.3 


6 


1.60 ±0.17 


18±2 


7 


>1000 


>1000 


8 


>1000 


>1000 


9 


17 ± 1 


310±28 


10 


>1000 


>1000 


11 


>1000 


>1000 


12 


34 ± 1 


120 ± 15 


13 


8.0 ± 1.0 


38 ±2 



The compounds 7, 8, 10, and 11 did not show any inhibitory activity {Ki > 1000 |xM). A moderate 
potency {Kj about 20-60 |xM) was observed in the case of rhizochalin (1) and its aglycone (2). 
Somewhat less active towards both T. californica and human al nAChRs were pibocin (3) and 
aaptamine (12). Crambescidin 359 (9) has shown a twenty-fold higher activity on muscle-type 
receptor; in other words, it has some subtype specificity as revealed in competition with [ I]-aBgt. 
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Makaluvamine G (6) has shown a similar ten-fold preference towards T. californica nAChR. The most 
potent competitor of [^^^I]-aBgt for the "classical" binding site on both tested receptors was 
makaluvamine C (5): it showed Ki = 1.30 ± 0.05 |a,M on muscle-type receptor and 3.6 ± 0.3 f^M on 
human al nAChR. 

Summarizing these results on the activity of the studied compounds from sponges and ascidians, we 
ascertain that there are substances with sound cholinergic properties in these marine sources acting on 
both muscle-type and a7 neuronal nAChRs. Their mode of action and affinities are close to those for a 
number of well-known cholinergic ligands from alkaloids (nicotine, lobeline, epibatidine, 
D-tubocurarine, DHpE) [41], to different a-conotoxins [42-44], some snake venoms components 
(waglerins, weak toxins) [45,46] and native modulators [47] showing similar micromolar potency to 
the muscle and a7 nAChR subtypes evaluated in competition with ['^^I]-aBgt. 

2.4. Electrophysiological Analysis of Effects on Functional Activity of Nicotinic Receptors 

To compare the affinity found in binding studies with possible functional activity, 
electrophysiological measurements were performed. We applied 10 |j,M solutions of tested compounds 
to murine muscle nAChR and to human al nAChR expressed in Xenopus laevis oocytes. Figure 5 
represents the inhibition activities of tested compounds in percents of maximal inhibition achieved by 
15 min application of 1 |^M a-bungarotoxin. 

Figure 5. Relative inhibition of agonist-evoked current by 10 |j.M compounds on murine 
muscle-type nAChR (a) and on human al nAChR (b). Numbering the compounds in x axis 
corresponds to the numbering in Figure 1. Asterisks indicate significant {p < 0.05, 
according to Student's test) differences between inhibition effects on murine muscle- and 
human al nAChRs for compounds (4), (6) and (11). 

a 



0 1 2 a 4 S 6 7 6 9 10 11 12 13 




0 1 2 3 4 5 6 7 8 9 10 11 12 13 



b 

0 1 2 3 4 5 6 7 8 9 10 11 12 13 

120 -j ' ' - — ' ■ ■ ■ ■ ' ■ ■ ' ■ ■ ■ — ^ 




0 1 2 3 4 5 6 7 8 9 10 11 12 13 
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Compounds 1, 2, 6, and 13 blocked muscle nAChR at given concentration almost completely, while 
compounds 3, 4, 9, and 12 showed moderate activity. A weaker but still reproducible effect was shown 
by compound 5, whereas only extremely weak, if any, activity was observed for compounds 7, 8, 10 
and 11. It is noticeable that the effects of compounds with a long non-polar carbon chain (1, 2 and 13) 
were very slowly reversible, despite their relatively good solubility in water as seen in Figure 6a for 
rhizochalin (1) (for compounds 2 and 13 data not shown). All the tested compounds, with the 
exception of makaluvamine G (6), varacin (4) and to a lesser extent compound (11), showed similar 
effects both on al and muscle nAChR at a given concentration. Although 10 |xM makaluvamine G (6) 
blocked muscle nAChR almost completely, like slowly-reversible rhizochalin and its aglycone (1, 2) 
and monanchocidin (13), its inhibition has a fast off-set rate and the receptor recovered in 5 min from 
the end of compound application (Figure 6b). Interestingly, 10 \iM varacin showed a strong inhibition 
of muscle nAChR, but not al nAChR. A difference in rates of inhibition of muscle and al nAChRs 
was observed also for compounds (6) and (11). However, even statistically significant differences in 
rates of inhibition should not be interpreted as a sign of selectivity towards one or another receptor 
subtype without further and more thorough studies. 

Figure 6. Electrophysiological measurements of 10 |xM rhizochalin (1) (a) and 
makaluvamine G (6) (b) activity on muscle nAChR. Black and grey rectangles represent 
application of acetylcholine and tested compound, respectively. From left to right: control 
acetylcholine-evoked current, acetylcholine-evoked current in the presence of tested 
compound and acetylcholine-evoked current after 15 min of wash out. Current inhibition 
dose-response curves (c): concenfrations of tested compounds from I to 100 |xM were used 
to evaluate IC50 (values placed on figure according to numbering from Figure 1). The 
sjmibols used for the respective compounds are the same as in Figures 2 and 4. 




Ig [ compound, M] 



We obtained IC50 values on muscle nAChR for some of the most potent compounds (Figure 6c). As 
can be seen, aglycone of rhizochalin (2) and makaluvamine G (6) are the most potent inhibitors of 
mouse muscle-type nAChR. The differences in the affinity and specificity for some compounds (1, 2, 3, 5) 
found in comparison of the data obtained by radioligand analysis and electrophysiology might depend 
on the intrinsic features of these two methods. For the former, a crucial factor is potent and virtually 
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irreversible binding of [ I]-aBgt to 2 binding sites on muscle-type nAChR and to 5 sites on a7 
nAChR. On the other hand, in electrophysiology experiments the choice and concentration of agonist 
are of importance. 

Summarizing the results of electrophysiology experiments, we can conclude that a number of 
analyzed compounds from marine sponges and ascidians possess functional (inhibitory) activities at 
micromolar concentrations. The majority of active compounds are not selective in respect to the 
muscle or neuronal a7 nAChRs, but some of them appear to be relatively selective towards studied 
receptor subtypes. However, fijrther experiments will reveal whether they might act on heteromeric 
neuronal nAChRs or other Cys-loop receptors. 

3. Experimental Section 

3.1. Isolation of Individual Compounds 

Purity of all the obtained samples were confirmed by TLC, HPLC, as well as by measurements of 
physical constants, while their structures were established mainly using NMR and MS data at 
comparison with those reported in literature [20-36]. 

3.2. Computer Modeling 

Avogadro software was used to generate energy minimized structures of ligands [48]. Molecular 
dockings were performed using Autodock 4.2/Adt software [49]. Keeping in mind that some recent 
studies [50] discussed spatial rearrangements in the Cys-loop proteins bound to diverse types of ligands, 
including partial and full agonists and antagonists, we decided to perform docking to those AChBP 
structures which were observed in complexes with HEPES, a-cobratoxin and dihydro-P-erythroidine 
(PDB 1UX2, 1YI5 and 4ALX, respectively), after removing these ligands. Initially, only those in silico 
structures involving the binding sites for agonists and competitive antagonists were taken into account. 
Search space was centered at the C-loop, x/y/z dimensions were 120/70/50 where x is coUinear with 
symmetry axis of the receptor and y is tangential to C-loop. The following Autodock Lamarckian 
genetic algorithm parameters were set: 200 dockings, evaluation number — long. Results were analyzed 
in respect of predicted affinity and compatibility with the aBgt binding site. LigPlot+ software [51] 
and PyMOL™ 1.2r2 (Molecular Graphics System, DeLano Scientific LLC, San Carlos, CA, USA) 
were used to analyze docked structures. 

For rhizochalin, aglycon of rhizochalin and monanchocidin (1, 2, 13, respectively) we decided to 
dock only part of the molecule to overcome Autodock' s limitations to the number of flexible bonds. 
For rhizochalin and its aglycon the part from the oxo-group to the closest sphingosine head was used. 
The pentacyclic (crambescidine-like) part of monanchocidin represented this ligand in the docking. 

3.3. Radioligand Assay 

\q competition experiments with [ I] -aBgt, all the compounds (in concentration range of 
0.03-1000 |xM) were pre-incubated 2-3 h at room temperature with the L. stagnalis AChBP (the final 
concenfration of 2.4 nM) in 50 [xL of buffer A (phosphate-buffered saline, 0.7 mg/mL of bovine serum 
albumin, 0.05% Tween 20, pH 7.5) or with the GH4C1 cells (6.5 \ig of total protein with final 
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concentration of 0.4 nM of toxin-binding sites) or Torpedo californica electric organ membranes (final 
concentration 1.25 nM of toxin-binding sites) in 50 |^L of buffer B (20 mM Tris-HCl buffer, 1 mg/mL 
of bovine serum albumin, pH 8.0). After that [ l]-aBgt was added to L. stagnalis AChBP, GH4C1 
cells or membranes to final concentration 0.1-0.2 nM and the mixtures were additionally incubated for 
5 min. Binding was stopped by rapid filtration on double DE-81 filters (Whatman, Maidstone, UK) 
pre-soaked in buffer A (for AChBP) or on GF/C filters (Whatman) pre-soaked in 0.25% 
polyethylenimine (for GH4C1 cells or membranes), unbound radioactivity being removed from the filters 
by washout (3x3 mL) with the buffers A and B, respectively. Non-specific binding was determined in 
all cases using 2-3 h pre-incubation with 10 |xM a-cobratoxin. 

The binding results were analyzed using ORIGIN 7.5 (OriginLab Corporation, Northampton, MA, 
USA) fitting to a one-site dose-response curve by Equation: 



where IC50 is the concentration at which 50% of the binding sites are inhibited and n is the Hill 
coefficient. The Ki values were calculated from the experimentally measured IC50S by the method of 
Cheng & Prusoff according to Equation: 



where L is concentration of free radioligand and Kd is dissociation constant for [ I]-aBgt to the 
respective target measured under similar assay conditions. Under radioligand assay conditions used the 
measured Kd values were 0.032 ± 0.003 nM for T. californica nAChR, 0.88 ± 0.19 nM for human a7 
nAChR and 3.5 ± 0.6 nM for L. stagnalis AChBP. 

3.4. Electrophysiology Measurements 

Oocytes were obtained from benzocaine anesthetized Xenopus by surgically removing a part of the 
ovarium. Single oocj^es were isolated and transferred to ND96 electrophysiology buffer. To achieve 
nAChR expression, each oocyte was injected with 1-5 ng of plasmid DNA containing the respective 
subunit genes (CHRNAl, CHRNBl, CHRND, CHRNE and CHRNA7) under control of CMV-promoter. 
Current recordings were performed after 36-48 h of incubation at 18 °C in ND96 solution. 

All recordings were performed using the Turbo TEC-03X amplifier (npi electronic GmbH, Tamm, 
Germany). 3M KCI-filled recording electrodes with a resistance of about 0.1-0.5 MQ were used. To 
achieve a minimum of compound consumption we used a hand-made recording chamber with 
approximately 25 |iL volume. Membrane potential was clamped at 50-70 mV. During the experiment 
each oocyte was perfused with ND96 solution until the resting current reached a steady state. 
Then short pulses (approximately 10 s) of agonist (50 |j,M nicotine in the case of a7 or 10 |xM 
acetylcholine in the case of muscle nAChR) were applied to the oocyte and then washed out with 
confrol buffer until the current returned to baseline. If the amplitude of agonist-evoked currents did not 
change after 5 min wash, oocj^es were considered suitable for fijrther experiments. Applications of 
agonist were performed every five minutes. Amplitudes of current responses to agonist in the presence 
of tested marine compounds were compared to the preceding agonist-induced current response. Each 
compound was tested on three different oocj^es expressing each type of receptor (muscle or a7). 
In preliminary experiments the marine compounds were tested at a concentration of 10 |xM, and 



% response = 100/{1 + ([toxinJ/ICso)"} 



(1) 



Ki = lCs(J{\+LIKd) 



(2) 
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QtiPlot software was used to display the results. For the most active compounds further measurements 
with concentrations of tested ligand in the range 1.25-20 \iM were done, and IC50 values were 
calculated using Origin 7.5 software (OriginLab, Northampton, MA, USA). 

4. Conclusions 

In this paper, we described the interaction between nicotinic acetylcholine receptors and several low 
molecular weight naturally-occurring compounds from marine sponges and ascidians. Molecular 
docking, radioligand binding and functional electrophysiological tests were done to determine potency 
of these products. Makaluvamines C (5) and G (6) were the most active in radioligand binding assay, 
and makaluvamine G (6) as well as crambescidine 359 (9) revealed most evident selectivity towards 
muscle-type nAChR. In electrophysiological tests rhizochalin (1), aglycone of rhizochalin (2), pibocin (3), 
varacin (4), makaluvamine G (6), crambescidine 359 (9) and monanchocidin (13) demonstrated 
highest potency (in the low micromolar range) on mouse muscle nAChR. The same compounds, with 
the exception of (4) and (6), showed the similar affinities on human a7 nAChR. Our results indicate 
that some of the earlier described pleiotropic activities of the analyzed marine products may be 
associated with their action on nicotinic acetylcholine receptors. On the other hand, the activities 
revealed here by radioligand analysis and electrophysiology, as well as computer modeling in 
comparison with the crystalline structures of complexes of nicotinic ligands, may be useful for 
designing drugs against diseases associated with malfunctioning of the cholinergic system (myasthenia 
gravis, psychiatric and neurodegenerative diseases). 
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